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HIGHLIGHTS 


•  Robust  condition  monitoring  of  Li-Ion  cell  using  multiple  model  adaptive  estimation. 

•  Equivalent  circuit  based  model  was  updated  with  a  nonlinear  OCV-SOC  relationship. 

•  The  OCV-SOC  equation  was  obtained  via  curve  fitting  of  the  experimental  data. 

•  The  model  bank  includes  a  normal  cell  and  two  distinctively  over-discharged  cells. 

•  RLS  method  was  used  to  identify  the  model  parameters  from  measured  data. 
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In  this  paper,  a  model  based  condition  monitoring  technique  is  developed  for  lithium-ion  battery  con¬ 
dition  monitoring.  Here  a  number  of  lithium-ion  batteries  are  cycled  using  two  separate  over  discharge 
test  regimes  and  the  resulting  shift  in  battery  parameters  is  recorded.  The  battery  models  are  constructed 
using  the  equivalent  circuit  methodology.  The  condition  monitoring  setup  consists  of  a  model  bank 
representing  the  different  degree  of  parameter  shift  due  to  overdischarge  in  the  lithium  ion  battery. 
Extended  Kalman  filters  (EKF)  are  used  to  maintain  increased  robustness  of  the  condition  monitoring 
setup  while  estimating  the  terminal  voltage  of  the  battery  cell.  The  information  carrying  residuals  are 
generated  and  evaluation  process  is  carried  out  in  real-time  using  multiple  model  adaptive  estimation 
(MMAE)  methodology.  The  condition  evaluation  function  is  used  to  generate  probabilities  that  indicate 
the  presence  of  a  particular  operational  condition.  Using  the  test  data,  it  is  shown  that  the  performance 
shift  in  lithium  ion  batteries  due  to  over  discharge  can  be  accurately  detected. 
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1.  Introduction 

Lithium  ion  (li-ion)  batteries  are  the  electrochemical  energy 
source  of  choice  today.  A  typical  li-ion  rechargeable  cell  with 
lithium  metal  oxide  based  positive  electrode,  graphitic  carbon 
negative  electrode,  and  lithium  conducting  organic  electrolyte,  of¬ 
fers  great  advantages  over  other  battery  chemistries  [1,2].  With 
major  advantages  of  high  energy  density,  safer  chemistries,  low 
self-discharge,  longer  cycle  life,  broad  temperature  application,  the 
li-ion  batteries,  availability  in  different  form  factors,  are  used  in  a 
range  of  applications  like  consumer  electronics,  automotive,  space 
exploration,  and  medical  implants  1,3  ,  to  name  a  few.  With  these 
growing  applications  in  mind  the  health  of  the  Li-ion  battery  be¬ 
comes  a  critical  factor  in  the  combined  system  functionality  of  the 
device. 
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Faults  occurring  in  the  li-ion  battery  can  be  attributed  to  number 
of  factors  either  individual,  like  structural  failure,  failure  of  thermal 
management,  or  more  likely,  a  combination  of  factors  involving 
manufacturing  defects,  over  charge,  over  discharge,  and  short  cir¬ 
cuit.  The  currently  available  li-ion  battery  safety  devices  can  be 
broadly  divided  into  internal  and  external  protection.  The  internal 
elements  are  implanted  on  the  battery  cell  and  provide  protection 
against  over  current,  high  temperature,  high  pressure,  over  charge, 
and  over  discharge.  The  commonly  used  internal  safety  elements 
are  the  polymeric  positive  temperature  coefficient  (PPTC),  charge 
interrupt  devices  (CID),  and  the  printed  circuit  boards  [4,5  .  The 
external  protection  elements  ensure  the  battery  safety  under  over 
charge,  overdischarge,  and  shorting  through  the  use  of  devices  like 
protection  diodes,  dedicated  battery  charging  integrated  circuit 
elements,  temperature  sensors,  and  more.  Some  of  these  devices 
are  resettable  while  others  are  one  time  use,  which  later  renders 
the  battery  useless.  Furthermore,  these  protective  devices  fail  to 
provide  the  user  with  any  information  regarding  the  condition  of 
the  battery,  extent  of  the  fault,  fault  identification,  and  battery 
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health  prognosis.  Critical  insight  into  the  battery  health  can  be 
obtained  by  implementing  analytical  redundancy.  It  involves 
reconstructing  the  process  behavior  on-line  by  using  models  which 
mimic  the  actual  process  under  study.  The  application  of  analytical 
redundancy  for  li-ion  battery  condition  monitoring,  fault  detection 
and  diagnosis  results  in  better  understanding  of  the  battery  dy¬ 
namics,  and  hence  contributes  towards  safer  batteries,  which  leads 
to  the  overall  system  safety. 

The  application  of  fault  detection  and  diagnosis  on  li-ion  bat¬ 
teries  is  not  new,  extensive  work  in  this  area  has  been  done  by 
researchers  with  focus  on  different  faults  and  related  techniques. 
The  work  on  li-ion  battery  fault  detection  is  based  primarily  on  the 
state  estimation,  empirical  techniques,  parameter  identification, 
data  driven  methods  and  others. 

Substantial  work  in  the  field  of  fault  detection  and  diagnosis, 
and  prognosis  in  li-ion  battery  using  data  driven  methods  has  been 
carried  out  by  Saha  et  al.  6-8].  Related  research  by  using  support 
vector  machine  algorithm  for  state  of  health  (SOH)  and  remaining 
useful  life  (RUL)  was  recently  carried  out  by  Nuhic  et  al.  [9]  and 
Wang  et  al.  10].  These  methods  involve  the  application  of  classi¬ 
fication  and  regression  algorithms  found  under  the  paradigm  of 
machine  learning.  In  Ref.  [11],  the  author  uses  AC  impedance 
spectroscopy  (IS)  along  with  auto  regressive  moving  average 
(ARMA),  neural  network,  and  fuzzy  logic  techniques  for  parameter 
identification,  estimation  and  eventually  battery  prognosis.  Data 
driven  techniques  do  not  require  in  depth  knowledge  of  the  battery 
and  its  underlying  mechanisms,  hence  their  implementation  does 
not  involve  expert  knowledge  of  the  process  under  study.  The 
biggest  hurdle  in  using  data  driven  methods  can  be  attributed  to 
the  computational  expensiveness  and  requirement  of  extensive 
data  for  training,  and  the  time  involved  in  learning. 

A  combination  of  rule  based  signal  monitoring  and  probability 
based  Li-ion  battery  fault  detection  and  diagnosis  was  explored  by 
Xiong  et  al.  [12],  these  methods  rely  heavily  on  the  thermal  sig¬ 
natures  of  the  battery  which  in  turn  depend  on  the  rate  of  charge/ 
discharge  applied  on  the  cell.  Further,  there  is  little  information 
regarding  the  initial  state  of  the  cell  under  test;  as  it  is  difficult  to 
achieve  an  over  discharge  cell  failure  in  LiFeP04  cell  chemistries 
after  two  cycles.  In  Ref.  [13  ,  the  open  circuit  voltage  (OCV)  is 
analyzed  along  with  model  based  approach  to  detect  the  cell 
nominal  capacity  fade  due  to  cycling.  This  technique  gives  good 
results  for  offline  applications  where  the  load  can  be  disconnected 
and  there  is  enough  time  to  accurately  access  the  OCV  of  a  given 
cell. 

Sate  estimation  involves  the  evaluation  of  the  state  of  the  bat¬ 
tery,  while  the  choice  of  technique  can  differ  based  on  the  re¬ 
quirements,  the  aim  is  to  access  the  information  related  to  the  Li- 
ion  battery  that  is  not  readily  available  through  measurement 
[14].  The  choice  of  state  variable  depends  on  the  model  of  the 
system,  but  for  Li-ion  batteries,  SOC  among  others  is  a  natural 
candidate.  Application  of  Luenberger  observers  (LO)  for  fault 
detection  and  diagnosis  can  be  found  in  Ref.  15],  here  the  authors 
implement  fault  diagnosis  on  a  string  of  Li-ion  batteries  using  a 
bank  of  reduced  order  observers.  LO  is  a  good  candidate  for  fault 
detection  and  diagnosis  in  systems  with  little  or  no  measurement 
noise,  but  with  presence  of  noise,  this  setup  will  face  inherent 
difficulties  especially  under  subtle  but  important  performance 
variation.  The  use  of  Kalman  filters  under  the  paradigm  of  observer 
based  fault  diagnosis  for  fault  detection  and  diagnosis  in  Li-ion 
batteries  is  given  in  Ref.  [16].  Where  the  optimal  filter  shows 
strong  robustness  to  noise  and  the  adaptive  nature  of  the  algorithm 
ensures  accurate  fault  detection. 

The  use  of  observer  based  fault  diagnosis  under  the  paradigm  of 
model  based  fault  diagnosis  offers  inherent  benefits  like  the 
decoupling  of  faults  of  interest  from  other  faults,  and  minimizing 


the  effects  of  unknown  disturbances  and  model  uncertainties  17]. 
These  advantages  are  further  utilized  in  the  multiple  model  adap¬ 
tive  estimation  (MMAE)  technique;  a  special  type  of  observer  based 
fault  diagnosis  technique.  MMAE  employs  a  Kalman  filter  bank  of  n 
filters,  where  one  observer  represents  the  healthy  condition  of  the 
process  being  monitored  while  the  remaining  n-1  observers 
represent  the  fault  conditions  of  the  process  [16,18,19  .  In  addition 
to  this  apparent  extension  to  the  single  observer  case,  MMAE  also 
provides  the  added  advantage  of  including  a  probabilistic  approach 
to  fault  detection  and  diagnosis. 

This  paper  is  organized  as  follows:  Section  2  describes  the  bat¬ 
tery  model,  Section  3  describes  the  model-based  fault  diagnosis 
using  nonlinear  observers  for  residual  generation  and  probability 
evaluation.  Section  4  discusses  the  design  of  the  experiment,  and 
Section  5  provides  the  discussion  of  the  results  obtained.  The 
conclusion  of  the  work  is  captured  in  Section  6. 

2.  Battery  modeling 

Li-ion  batteries  can  be  modeled  using  different  techniques 
namely  electro  chemical,  neural  networks,  empirical,  experimental 
and  equivalent  circuit  [20,21  .  The  choice  of  modeling  technique  is 
a  tradeoff  between  capturing  cell  dynamics  and  computational 
demand.  For  real  time  application  the  equivalent  circuit  model 
approach  is  adopted  because  it  gives  good  representation  of  cell 
dynamics  while  maintaining  low  computational  resource  usage. 

The  Li-ion  battery  can  be  modeled  as  a  third  order  system  using 
lumped  electrical  elements  like  resistors  and  capacitors.  The 
equivalent  circuit  model  is  shown  in  Fig.  l.where,  Rb  is  the  ohmic 
resistance,  which  accounts  for  the  limited  conductance  of  the 
metallic  contacts,  inter  cell  connections,  electrode  material  and  the 
bulk  electrolytic  resistance  to  electron  and  ion  migration  [1,22], 
constant  phase  element  (CPE)  C  and  resistance  R  are  used  to  model 
the  distribution  of  reactivity  depicting  the  local  property  of  the 
electrode.,  charge  transfer  resistance  Rct  and  double  layer  capaci¬ 
tance  Cdi  represent  the  interfacial  impedance  of  the  cell  [23]  and 
Vocv  represents  the  battery  cell  OCV.  The  CPE  captures  the  distri¬ 
bution  of  reactivity  at  the  electrodes  which  can  be  attributed  to 
variation  in  surface  properties.  The  impedance  function  of  the 
combined  RC  pair  is  given  by  Refs.  [23,24  , 

zcpeM  =  1  —  (joj)aQR  (1) 

where,  a  is  the  depression  factor  associated  with  the  CPE  and  is 
assumed  to  be  unity.  As  a  result  Qcan  be  replaced  by  C  and  the  CPE 
then  behaves  like  a  normal  capacitor  22,23]. 

The  circuit  parameters  depend  on  the  SOC,  temperature  and 
capacity  fade  effects  [25  .  For  this  study,  parameter  dependency  on 
these  factors  is  assumed  to  be  small.  The  effect  of  non-linear 
element  in  the  equivalent  circuit  namely  Warburg  impedance 
representing  the  diffusion  phenomenon  is  considered  to  be  negli¬ 
gible  [26  . 
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Fig.  1.  Li-ion  battery  equivalent  circuit  model. 
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The  function  /( SOC)  represents  the  non-linear  function  that 
maps  the  relationship  between  the  OCV  and  the  SOC  of  the  li-ion 
battery.  This  relationship  is  given  by  the  classical  OCV-SOC  curve, 
as  shown  in  Fig.  2. 

The  OCV-SOC  data  was  recorded  from  a  sample  LiFeP04  battery 
cell  tested  at  room  temperature  at  the  Energy  Systems  and  Power 
Electronics  Laboratory  (ESPEL)  at  IUPUI.  Interested  readers  are 
directed  to  Ref.  27]  for  further  reading  on  OCV-SOC  relationship 
evaluation.  Unlike  the  almost  linear  trend  shown  by  lead  acid 
batteries  [1  ,  the  OCV-SOC  profile  for  li-ion  battery  shows  non¬ 
linear  behavior  with  a  relatively  flat  trend  between  20  and  80% 
SOC.  In  cases  where  the  OCV  is  used  as  an  indicator  of  battery  SOC,  a 
small  error  in  OCV  evaluation  can  result  in  appreciable  error  in  the 
resulting  SOC.  In  addition,  temperature  along  with  age  also  has  an 
effect  on  the  OCV-SOC  trend,  and  can  cause  it  to  shift. 

The  OCV-SOC  relationship  in  Fig.  2  can  be  captured  by  a  poly¬ 
nomial  with  varying  degree  and  coefficients,  resulting  in  different 
degrees  of  fit.  A  ninth  degree  polynomial  is  found  to  give  the  best  fit 
and  is  selected.  The  non-linear  function  /(SOC)  is  then  given  by, 

/(SOC)  =  cig(SOC)9  +  a8(SOC)8  +  a7(SOC)7  +  a6(SOC)6 

+  a5(SOC)5  +  a4(SOC)4  +  a3(SOC)3  +  a2(SOC)2 
-T  Oi  (SOC)  +  Qq  (2) 


where  VCdl  is  the  voltage  across  the  double  layer  capacitor  Cdi. 

The  SOC  is  defined  as  the  ratio  of  the  remaining  capacity  to  the 
fully  charged  nominal  capacity  of  the  battery  [28],  and  is  given  by, 

t 

SOC (t)  =  SOC(O)  +  [  (5) 

J  En 

0 

where  SOC(O)  represents  the  initial  state  of  charge,  Cn  represents 
the  battery  cell  nominal  capacity  in  Ampere  hour,  ri  is  the  coulomb 
efficiency  given  by, 

f  1 ,  charging 
1  \  0.98,  discharging 

In  discrete  time,  Eq.  (5)  can  be  given  by, 

SO C(k)  =  SO C(k  -  1 )  +  vAlrL(k>  (6) 

where  k  is  the  time  variable,  and  A k{k)  =  Jj,(k)-t(k-l). 

The  discrete  time  versions  of  Eqs.  (3)  and  (4)  can  be  obtained  by 
using  zero-order  hold  (ZOH)  process  [29  .  The  transformed  equa¬ 
tions  are  given  by, 


where  ag  =  0.0385,  ag  =  -0.01936,  a-j  =  -0.169,  as  =  0.06142, 
<35  =  0.2328,  <34  =  -0.05715,  <33  =  -0.08321,  <32  =  0.0005257, 
eq  =  0.03205,  a0  =  3.297. 

From  the  equivalent  circuit  in  Fig.  1,  using  Kirchhoff's  voltage 
law,  the  rate  of  change  in  voltage  across  the  capacitor  C  is  given  by, 


Vc^k 


RC  C 


where  Vc  is  the  voltage  across  the  capacitor  C,  and  IL  is  the  battery 
load/charge  current.  The  sign  convention  used  in  this  study  con¬ 
siders  the  negative  sign  of  k  as  discharging  while  positive  sign  as 
charging.  The  rate  of  change  in  voltage  across  the  double  layer 
capacitor  is  given  by, 
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Fig.  2.  Experimental  OCV-SOC  curve  for  LiFeP04  battery  cell. 
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The  battery  terminal  voltage  at  any  time  sample  k,  is  given  by, 


Vt(fc)  =/( SOC)  -  IL(k)Rb  -  Vc(k)  -  VCdl(k)  (9) 

where  Vt  is  the  battery  terminal  voltage,  and  kRb  is  the  voltage  drop 
across  the  bulk  resistance  % 

The  state  vector  is  given  by  xk  =  [  SOC (k)  Vc  (/<)  (/<)  ] T  and 

the  state  equation  for  a  non-linear  time  invariant  system  in  discrete 
time  is  given  by, 


X(fc)  =  g(*/<-l,  Uk-P+Wfc-!  (10) 

and  the  output  equation  is  given  by, 


z(k)  =  h(xk,  uk)  +  vk  (11) 

where  g  and  h  are  continuously  differentiable  non-linear  functions, 
w  is  the  process  noise  with  zero  mean  and  variance  of, 

£{w„[/K[m]}  =  {^’  l=%  (12) 

and  v  is  the  measurement  noise,  independent  from  w,  with  zero 
mean  value  as, 

E{v„[l]vTn[m]}  =  {J  l  =  ™  (13) 

Q.  and  R  are  the  process  and  measurement  noise  variances 
respectively.  The  process  and  measurement  white  Gaussian  noise  is 
generated  using  the  polar  method  [30  . 
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From  Eqs.  (6)— (11),  the  functions  g  and  h  are  given  by, 


S(fe-l)  = 


SOC(k-  1)  + 


IL(k  —  1) 


Vc{k—  1)  +  J? 


I(k-  1) 


e  fictQi 

VCd](k  -  1)  +  Rct 

1  - 

Q  ftctQl 
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/ 

(14) 


and 

m  =/( SOC)  -  IL(k)Rb  -  Vc(k)  -  VCdl(k)  (15) 

3.  Model  based  fault  diagnosis 

The  choice  of  diagnosis  method  is  based  on  a  combination  of 
factors  like  cost,  range  of  operation  and  complexity  of  process  [17]. 
Software/analytical  redundancy  fault  diagnosis  is  particularly 
important  because  it  can  be  implemented  at  a  low  cost,  incorporate 
process  complexity  and  encompass  range  of  system  operation. 
Model  based  fault  diagnosis  is  an  integral  part  of  software  redun¬ 
dancy  based  fault  diagnoses  and  is  based  on  the  process/physical 
model  of  the  system.  In  the  model  based  fault  diagnosis  scheme, 
the  process  model  is  the  mathematical  representation  of  the  sys¬ 
tem.  While  this  model  should  be  able  to  show  good  agreement  with 
the  system  dynamics  of  the  process  it  should  more  importantly 
conform  to  the  fault  carrying  dynamics  of  the  process.  Also  for  real 
time  applications,  a  computationally  inexpensive  system  repre¬ 
sentation  will  be  advantageous  as  computational  resources  are 
limited;  this  consideration  takes  all  the  more  precedence  in  mobile 
application.  The  measured  output  from  the  process  is  compared 
with  the  model  response  to  generate  the  residuals.  The  perfor¬ 
mance  of  the  model  based  fault  diagnosis  is  further  improved  by 
replacing  the  process  model  and  subsequent  residual  generation 
process  with  state  observers  17,31,32  ,  hence  called  the  observer 
based  fault  diagnosis.  This  fundamental  shift  in  residual  generation 
process  enables  the  removal  of  errors  in  the  residual  signal  due  to 
faulty  initial  condition,  and  unknown  disturbances. 

The  general  configuration  of  model  based  fault  diagnosis  with 
Kalman  filter  state  observer  is  as  shown  in  Fig.  3.  The  optimal  filter 


works  in  parallel  with  the  battery,  and  based  on  the  measurements 
of  load/charge  current  and  the  terminal  voltage,  it  estimates  the 
states  of  the  battery  model.  The  states  are  further  utilized  to 
generate  the  estimated  terminal  voltage  Vt.  The  fault/condition 
information  carrying  residuals  are  generated  when  the  estimated 
terminal  voltage  is  compared  with  the  measured  terminal  voltage 
[16]. 

In  the  case  of  non-linear  process  models  the  Kalman  filter  in  the 
observer  based  fault  diagnosis  is  replaced  by  the  extended  Kalman 
filter. 

3  A.  Extended  Kalman  filter  design 

The  Extended  Kalman  filters  are  used  for  state  estimation  of 
non-linear  systems  by  linearizing  around  the  current  mean  and 
covariance.  When  applied  to  the  non-linear  system  of  (10),  (11), 
(14),  and  (15),  the  time  update  equations  are  given  by  Refs.  [33,34], 

Pk  =  CkPk-,CTk  +  Qjc_i  1  ; 

and  the  measurement  update  equations  are  given  by, 

Kk  =  PkHTk(HkPkHTk+Rky 1 

xk  =  xk  +Kk(yk-h(xky  07) 

Pk  =  (/  -  KkHk)Pk 

where  xk  represents  predicted  state  based  on  the  function  g  eval¬ 
uated  at  the  previously  estimated  state  and  available  input,  Py  is 
the  a  priori  estimate  error  covariance,  K/<  is  the  Kalman  gain,  xk  is 


Fig.  3.  Observer  based  fault  diagnosis. 
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the  updated  state  estimate  and  P/<  is  the  updated  covariance  esti¬ 
mate.  The  state  transition  matrix  G/<  and  observation  matrix  H/<  are 
evaluated  at  each  step  and  are  given  by, 


*k- 1 5 Uk- 1 


*k  > 


The  estimated  terminal  voltage  is  given  by, 


(18) 


Vt{k)  =  h(xk,uk)  (19) 

As  mentioned  earlier,  the  residual  signal  is  obtained  by  sub¬ 
tracting  the  estimated  terminal  voltage  signal  from  the  measured 
terminal  voltage  of  the  battery  and  is  given  by, 


r(k)  =  Vt(k)  -  Vt(k) 


(20) 


where  r  is  the  fault/condition  residual,  Vt  is  the  measured  terminal 
voltage. 

The  single  observer  case,  as  shown  in  Fig.  3  can  be  extended  to 
multiple  observers  with  added  advantages.  Using  a  bank  of 
extended  Kalman  filters,  where  each  observer  represents  a  partic¬ 
ular  fault/operational  condition  of  the  battery  allows  for  robust  li- 
ion  battery  monitoring  covering  a  wide  variety  of  operational 
conditions. 


3.2.  Multiple  model  adaptive  estimation 


monitoring.  The  general  layout  of  MMAE  when  applied  to  non¬ 
linear  battery  model  with  n- 1  distinct  operational  conditions  is 
as  shown  in  Fig.  4. 

The  residual  generation  process  involves  non-linear  observers 
which  run  parallel  to  the  battery  with  all  the  n  filters  receiving 
the  same  input  of  load/charge  current  k  and  terminal  voltage 
measurement  Vt.  The  battery  condition  is  continuously  moni¬ 
tored  and  the  n  residuals  are  updated  at  each  time  instant.  The 
post  processing  part  involving  the  residual  processing,  decision 
logic  and  subsequent  fault/condition  declaration  takes  place  at 
the  conditional  probability  density  evaluator  block.  At  the  con¬ 
ditional  probability  evaluator,  each  of  the  n  fault/operational 
hypothesis  is  assigned  a  probability  of  being  true.  The  probabil¬ 
ities  range  between  0  and  1  and  the  sum  of  all  the  n  probabilities 
is  equal  to  1. 

The  earliest  use  of  conditional  probability  density  evaluation 
and  weighting  of  coefficients  can  be  found  in  Refs.  35-37],  where 
the  primary  emphasis  is  the  optimal  state  estimation  and  its  ap¬ 
plications.  With  respect  to  fault  detection  and  diagnosis  based 
implementation  of  Kalman  filter  and  conditional  probability  den¬ 
sity  evaluation,  extensive  work  has  been  done  by  Maybeck  et  al. 
[38-42]  and  Athans  et  al.  [43]  on  the  fault  detection  and  diagnosis 
of  aircraft  systems. 

The  conditional  probabilities  require  a  priori  samples  to 
compute  the  current  values  and  are  normalized  over  a  complete 
sum  of  conditional  probabilities  of  all  systems.  The  largest  condi¬ 
tional  probability  among  all  can  be  used  as  an  indicator  of  the  fault/ 
operational  hypothesis  being  true.  The  probability  for  the  nth  model 
at  time  sample  k  is  given  by, 


MMAE  is  a  special  type  of  observer  based  fault  diagnosis  tech¬ 
nique  as  it  employs  a  Kalman  filter  bank  (KFB)  of  n  filters,  where 
one  observer  represents  the  healthy  condition  of  the  process  being 
monitored  while  the  remaining  n-1  observers  represent  the  fault 
conditions  of  the  process  [16,19].  In  addition  to  this  apparent 
extension  to  the  single  observer  case,  MMAE  also  provides  the 
added  advantage  of  including  a  probabilistic  approach  to  condition 


fz(k)\a,Z(k-t)(^k  l®n>  ^k-\ )Pn(k  1 ) 
^2j=\fz(k)\a,Z(k-l)  (?k I®/?  ^k-\ )Pj( k  ~  1 ) 


(21) 


wher e/Z(/<)|aZ(/<_1)(Z/<|an,Z/<_1)  is  the  conditional  probability  density 
function  of  the  nth  model  considering  the  history  of  measurements 

Z(ti-i)  =  (zr(t1)...zr(ti_1)). 


Fig.  4.  MMAE  residual  generation  and  probability  evaluation. 
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The  conditional  probability  function  is  given  by, 

fz(k)\a,Z(k-l)  1 )  =  (3n  exp(°) 


Vt(k)=a9( SOC)9  + . . .  +a0(SOC)°  -IL(k)R„  -  Vc(k)  -  VQl (k)  (27) 

separating  the  unknown  parameters  from  the  known  signals,  the 
parametric  form  for  the  battery  model  is  given  byt, 


/  \ 

At 

vt(k)  = 

a9  ...  a0  Rt,  e~&  Rjl-e-^ 

e  RctCdi  Rct 
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_  e  PctQi 

V  ) 
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-i  T 
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(SOC(k))y  ...  (SOC(k))°  -IL(k)  -Vc(k- 1)  -IL(k- 1)  -vc  (k-1)  -IL(k - 1) 


dl 


(28) 


where 


1 


(2irf2  fn(k) 


1/2 


(23) 


l  is  the  measurement  dimension  and  equal  to  1,  since  VfeR. 


(24) 


where  rnjc  is  the  residual  signal  for  the  nth  model  at  time  sample  k. 
The  covariance  of  the  residual  signal  evaluated  at  each  sample  is 
given  by, 


^n,k^n,k\k^nj<  +  ^ 


(25) 


where  Cn,k  is  the  output  vector  for  the  nth  system  at  any  time 


sample  k.  For  the  non-linear  battery  model  Cn  k  =  ^ 


,  linearized 


%k\k 


output  vector  evaluated  at  the  current  estimated  state.  Pn  k\k  is  the 
updated  covariance  update  for  nth  model  evaluated  at  /cth  sample, 
R  is  the  measurement  noise  variance. 


4.  Design  of  experiments 

In  cases  where  the  battery  is  subjected  to  cycles  involving 
abusive  charge  or  discharge,  the  battery  undergoes  a  permanent 
change  in  its  performance.  The  shift  in  battery  performance  induces 
a  marked  variation  in  model  circuit  parameters  namely,  the  bulk 
resistance,  constant  phase  elements,  charge  transfer  resistance  and 
the  double  layer  capacitance,  from  their  healthy  battery  counter¬ 
parts.  The  parameter  values  of  the  battery  under  different  condition 
is  extracted  using  the  recursive  least  squares  technique  [44  . 

4.1.  Recursive  least  squares 

The  recursive  least  squares  (RLS)  system  identification  tech¬ 
nique  aims  at  fitting  battery  mathematical  model  to  a  sequence  of 
observed  battery  current  and  voltage  data  by  minimizing  the  sum 
of  the  squares  of  the  difference  between  the  observed  and 
computed  data  recursively  [44  .  The  system  being  identified  is 
represented  in  the  discrete  time  parameter  form  given  by, 

z(k)  =  e*T<p  (26) 

where  z  is  the  terminal  voltage  output,  0*  is  the  linear  vector  of 
unknown  parameters  that  is  being  identified,  cf>  is  the  vector  of 
previous  current  and  voltage  measurements. 

From  Eqs.  (9)  and  (2),  the  terminal  voltage  output  is  given  by, 


For  further  simplification,  some  of  the  parameters  can  be  lum¬ 
ped  together  as  shown  below, 


/ 


A^=e-&^=R[  1  -e~&  \  ,A2  =  e  R^1B2=Ra 
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At 


\~e  RctCdi 


/ 


(29) 


The  estimation  equation  is  given  by, 


d(k)  =  8(k  -  \)  +  P(k)<p(k)e(k) 


(30) 


where  P  is  the  covariance  matrix,  e  is  the  normalized  estimation 
error,  0  is  the  same  as  given  in  Eq.  (26).  At  k  =  0,  0( 0)  is  the  best 
initial  guess  on  the  parameters  to  be  identified. 

The  normalized  estimation  error  e  is  given  by, 


Vt(k)  -  djkfcpjk ) 
m2(k) 

where  m2  is  the  normalizing  signal  given  by, 


(31) 


m2(k)  =  cc  +  cpT(k)(p(k) 


(32) 


with  cc  >  0,  and  cceR 

The  covariance  matrix  P  is  updated  recursively  by  using  the 
following  equation, 


P(/<)  =  P(/<  -  1) 


P(k-  -l)<p(k)tp(k)T P(k  -  1) 
m2(k)  +  <p(k)TP(k  -  \)<p(k) 


(33) 


with  P( 0)  =P0  =  Pl>  0 


4.2.  Test  setup 


The  li-ion  battery  selected  for  this  study  was  A123  18650 
LiFeP04  (APR18650M  1A  3.3  V  1000  mAh)  from  A123  Systems 
(Cambridge,  MA)  [45].  The  battery  under  test  was  subjected  to 
nominal  charge/overdischarge  in  a  cyclic  fashion.  The  over¬ 
discharge  regime  is  based  on  Navy  overdischarge  cycle  [46]  and  a 
24  h  overdischarge  cycle.  In  the  Navy  over  discharge  cycle,  the  Li- 
ion  battery  is  discharged  at  maximum  suitable  discharge  rate  for 
25%  over  discharge.  The  charging  of  the  battery  is  carried  out  using 
a  standard  non-abusive  charge  regime.  The  battery  is  cycled  25 
times  using  this  discharge-charge  regime  and  the  critical  battery 
parameters  are  continuously  monitored.  In  the  24  h  over  discharge 
test  regime,  the  battery  is  discharged  at  a  suitable  discharge  rate 
until  the  SOC  reaches  zero.  To  maintain  near  zero  terminal  voltage, 
a  resistor  is  connected  across  the  terminals  for  the  duration  of  24  h. 
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Fig.  5.  Li-ion  battery  test  setup. 

The  charging  of  the  battery  is  then  carried  out  using  a  standard 
non-abusive  charging  profile.  The  24  h  over  discharge  test  cycle  is 
repeated  twice  and  the  battery  parameters  are  continuously 
monitored.  The  charging-discharging  and  data  acquisition  setup  is 
as  shown  in  Fig.  5. 

In  the  test  setup,  a  Mastech  programmable  DC  power  supply  is 
used  for  charging  of  the  battery,  powering  the  charge-discharge 
switching  circuit  and  the  battery  temperature  measurement  circuit. 
The  BK  Precision  programmable  DC  electronic  load  is  used  to 
implement  the  dynamic  discharge  of  the  battery  at  the  desired 
rates.  The  data  acquisition  is  carried  out  using  the  dSPACE  DS1103 
board  along  with  ControlDesk  and  the  data  post  processing  is 
implemented  in  MATLAB  R2009  environment. 

The  recorded  battery  current  and  voltage  data  at  the  end  of  both 
test  regimes  is  as  shown  in  Fig.  6.  The  signals  are  filtered  for  noise 
using  cascading  low  pass  filters  [47  . 

From  the  recorded  current  and  voltage  data  of  Fig.  6,  the 
parameter  values  for  the  battery  cell  before  and  after  the  nominal 


Table  1 

Identified  li-ion  battery  parameter  values. 


Parameter 

New  battery 

Navy  over  discharge 
cycled  battery 

24  h  Over  discharge 
cycled  battery 

a9 

0.49161 

0.748268 

-0.02105 

a8 

-0.08717 

-0.07862 

-0.04814 

a7 

-0.32599 

-0.42834 

-0.00855 

a6 

-0.27117 

-0.37816 

0.057027 

a5 

-0.03456 

-0.08725 

0.092889 

a4 

0.198062 

0.201331 

0.042196 

a3 

0.206745 

0.219284 

-0.10774 

a2 

-0.07741 

-0.10596 

-0.20669 

a1 

-0.00024 

-0.03705 

0.341718 

a° 

3.3558 

3.490569 

3.223314 

Rb 

-0.06492 

-0.05705 

-0.09248 

Ai 

0.241634 

0.607614 

0.216303 

Bi 

-0.06881 

-0.06389 

-0.09306 

A2 

0.241634 

0.034997 

0.216303 

b2 

-0.06881 

-0.06389 

-0.09306 

charge/over  discharge  cycles  are  estimated.  The  identified  system 
parameters  are  given  in  Fable  1. 

Using  the  identified  parameters,  multiple  battery  condition 
representing  models  can  be  formulated.  For  the  purpose  of 
validation,  the  load  current  applied  to  condition  monitoring 
setup  is  based  on  the  urban  dynamometer  driving  schedule 
(UDDS)  drive  cycle  profile  which  has  been  appropriately  scaled 
to  match  the  nominal  capacity  of  one  cell.  The  UDDS  drive  cycle 
is  accessed  from  Autonomie  [48]  and  is  represented  in  Fig.  7. 

Three  UDDS  cycles  are  run  back  to  back  with  the  total  testing 
duration  of  213  s.  While  the  load  current  profile  simulates  the 
actual  working  condition  of  the  system,  the  resulting  fault  proba¬ 
bilities  depend  more  on  the  zero  average  residual  signal  rather  than 
the  magnitude  of  the  load  current  16]. 

5.  Diagnosis  performance 

To  test  the  performance  of  the  condition  monitoring  setup,  a 
scenario  is  created  with  consecutive  change  in  the  battery  condi¬ 
tion.  The  total  simulation  time  of  213  s  is  divided  into  four  equal 


Load/charge  current 


0) 


O 


Terminal  voltage 


Fig.  6.  Current  and  voltage  profiles  for  healthy  and  nominal  charge/overdischarge  battery. 
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Fig.  7.  Battery  cell  load  current  profile. 


parts,  depicting  the  healthy  and  the  deteriorated  performance  of 
the  battery  due  to  overdischarge.  The  condition  change  scenario  is 
as  given  below, 

•  Zero  to  53.25  s:  healthy  battery  operation 

•  53.26-106.5  s:  battery  operation  after  two  24  h  overdischarge 
cycles 

•  106.51-159.75  s:  battery  operation  after  twenty  five  Navy 
overdischarge  cycles 

•  159.76-213  s:  healthy  battery  operation 

Once  the  operational  condition  is  diagnosed  correctly,  this  setup 
helps  to  check  the  effectiveness  of  the  condition  monitoring  algo¬ 
rithm  to  de-latch  itself  from  its  earlier  diagnosis  [49].  It  is  also 
assumed  that  only  one  type  of  fault  can  occur  in  the  system  at  any 
given  point  in  time.  The  starting  SOC  of  the  battery  is  considered  to 
be  at  70%  and  the  polarization  voltages  are  considered  to  be  zero. 


Simulated  &  estimated  terminal  voltage 


Simulated  &  estimated  terminal  voltage 


Fig.  9.  Simulated  and  estimated  terminal  voltages:  zero  to  53.25  s. 


The  UDDS  load/charge  current  and  the  terminal  voltage  mea¬ 
surement  embedded  with  the  changing  battery  condition  accord¬ 
ing  to  the  scenario  constructed  above  is  made  available  to  the  three 
non-linear  observers.  The  simulated  terminal  voltage  measurement 
and  observer  outputs  are  as  shown  in  Fig.  8. 

Fig.  9  below  shows  the  first  section  of  the  total  simulated  and 
estimated  terminal  voltages  in  Fig.  8. 

The  terminal  voltage  measurement  ym  depicts  the  healthy  bat¬ 
tery  operation  and  is  shown  by  the  continuous  blue  line  in  Fig.  9. 
For  the  duration  of  the  first  and  last  53.25  s,  the  estimated  terminal 
voltage  for  the  healthy  battery  system,  given  by  yNew  shows  good 
match  with  ym  while  the  estimated  terminal  voltages  for  the 
24  h  overdischarge  y24  hr od  -  and  Navy  overdischarge  yNavyoo  show 
marked  deviation. 

The  second  section  capturing  the  terminal  voltage  signals  from 
53.26  to  106.5  s  is  as  shown  in  Fig.  10. 

In  Fig.  10,  y24  hr0D  shows  good  match  with  ym  while  TNavyOD 
shows  large  deviation.  yNew  shows  large  deviation  from  ym  for 


Simulated  &  estimated  terminal  voltage 


Fig.  8.  Simulated  and  estimated  terminal  voltages. 


Fig.  10.  Simulated  and  estimated  terminal  voltages:  53.26-106.5  s. 
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Simulated  &  estimated  terminal  voltage 


Fig.  11.  Simulated  and  estimated  terminal  voltages:  106.51-159.75  s. 

higher  magnitudes  of  load/charge  current.  At  smaller  magnitudes 
of  load/charge  current,  yNew  is  close  to  ym  but  not  as  close  as 

y  24  hrOD- 

The  third  section  of  simulated  and  estimated  terminal  voltages 
is  as  shown  in  Fig.  11. 

In  Fig.  11,  y Navy od  shows  good  match  with  ym  while  the  other 
estimated  terminal  voltages  ^New  and  T24hrOD  show  large  de¬ 
viations  from  the  simulated  terminal  voltage  measurement. 

The  expected  behavior  shown  by  the  estimated  terminal  voltage 
signals  in  Figs.  9-11  leads  to  system  residuals  as  shown  in  Fig.  12. 
The  residuals  are  evaluated  using  Eq.  (20)  and  are  updated  at  each 
time  sample. 

During  the  first  and  the  last  53.25  s,  the  battery  heath  residual 
rnew  shows  a  zero  average  behavior  while  24  h  overdischarge  sys¬ 
tem  residual  r24  hr  od  and  Navy  overdischarge  system  residual  rNavy 
od  show  dynamic  behavior  within  the  same  time  durations.  This 
response  from  the  residuals  shows  that  the  healthy  system  opera¬ 
tion  is  accurately  captured.  Similarly,  from  53.26  to  106.5  s  T2a  hr  od 
shows  zero  average  behavior  while  other  residuals  show  dynamic 


- 1 - 1 - 

rs*yOO 

Time  [x  0.001]  (Sec)  ^5 


Conditional  probability  densities 


Time  [x  0  001]  (Sec)  x  ^q5 


Fig.  13.  Conditional  probability  densities  for  new  and  over  discharged  battery. 

behavior,  thus  confirming  the  presence  of  24  h  overdischarge 
condition.  With  accurate  state  estimations,  the  same  behavior  can 
be  seen  for  the  residual  rNavy  od- 

The  residuals  are  evaluated  using  the  probability  density  func¬ 
tion  and  the  normalized  probabilities  are  updated  at  each  time 
sample.  The  resulting  battery  condition  probabilities  are  as  shown 
in  Fig.  13. 

The  24  h  overdischarge  condition  was  injected  at  53.25  s,  which 
causes  the  battery  health  probability  Pnew  to  transition  from  1  to  0, 
indicating  the  non-existence  of  healthy  battery  condition.  The 
particular  battery  condition  is  further  indicated  by  the  probability 
^24hr  od  increasing  from  0  to  1.  At  the  same  time  a  spike  in  Navy 
overdischarge  probability  P^vy  od  can  be  seen,  this  is  because  the 
probabilities  can  take  some  time  to  transition  from  one  extreme 
state  to  another.  The  Navy  overdischarge  condition  is  correctly 
indicated  at  106.51  s,  where  the  condition  probability  PNavy  od 
transitions  from  0  to  1  and  P24hr  od  drops  from  1  to  0  while  Pnew 
remains  at  0.  Finally,  at  159.76  s,  the  healthy  battery  condition  is 
indicated  by  the  probability  Pnew  when  it  transitions  from  0  to  1.  At 


SOC  vanation 


Fig.  12.  System  residuals  for  healthy  and  over  discharged  battery. 


Fig.  14.  SOC  variation  for  healthy  and  over  discharged  battery  conditions. 
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the  same  time  PNavy  od  drops  from  1  to  0.  The  battery  conditional 
probabilities  show  good  and  accurate  behavior  by  detecting  the 
correct  battery  condition  at  the  correct  time  and  also  by  avoiding 
intermediate  loss  in  probability. 

The  SOC  is  one  of  the  elements  of  the  state  of  the  non-linear 
battery  model  and  contributes  greatly  towards  accurate  OCV  esti¬ 
mation  and  thus  the  terminal  voltage  estimation.  The  variation  of 
the  simulated  SOC  measurement  SOCm  using  coulomb  counting 
and  the  estimated  SOC  signals  is  as  shown  in  Fig.  14. 

For  accurate  condition  monitoring  of  the  li-ion  battery  it  is 
imperative  that  the  OCV  is  estimated  accurately  and  is  within 
reasonable  limits.  Inaccurate  OCV  estimation  results  in  loss  of 
critical  battery  condition  information  carried  by  the  polarization 
voltages.  To  better  approximate  the  actual  operation,  the  SOC  is 
bounded  by  0%  at  bottom  and  100%  at  the  top. 

In  Fig.  14,  during  the  first  and  the  last  53.25  s,  the  estimated 
healthy  battery  SOC  given  by  SOCnew  matches  with  SOCm  while  the 
estimated  24  h  overdischarge  SOC  given  by  SOC24hr  od  and  the  Navy 
overdischarge  SOC  given  by  SOCNavy  od  show  large  deviation.  From 
53.26  to  106.5  s,  SOC24hr  od  shows  good  match  with  SOCm  while 
SOCNavy  od  and  SOCnew  show  large  difference.  Finally  from  106.51  to 
159.75  s  SOCNavy  od  closely  follows  SOCm  while  SOCnew  and  SOC24hr 
od  show  apparent  deviation.  Using  the  extended  Kalman  filters,  the 
SOC  is  estimated  with  less  than  1%  accuracy  within  the  appropriate 
region  of  battery  condition. 

6.  Conclusion 

In  this  paper  an  observer-based  condition  monitoring  technique 
for  a  li-ion  battery  is  developed  and  validated.  Overdischarge  in  li- 
ion  battery  causes  change  in  the  operational  characteristics  of  the 
battery  which  can  be  identified  in  real  time  using  the  proposed 
condition  monitoring  technique.  The  non-linear  battery  model 
parameters  are  extracted  from  the  experimental  voltage  and  cur¬ 
rent  data  pertaining  to  healthy  and  over  discharged  batteries.  Bank 
of  extended  Kalman  filters  are  used  to  access  the  internal  dynamics 
of  the  battery  cell  and  to  estimate  the  internal  condition  of  the 
electrochemical  system.  The  effectiveness  of  the  condition  moni¬ 
toring  algorithm  was  tested  by  using  two  different  overdischarge 
regimes  on  the  li-ion  battery,  inducing  different  degree  of  battery 
condition  deterioration.  The  developed  condition  monitoring 
technique  showed  robust  performance  by  detecting  the  different 
degrees  of  overdischarge  on  the  li-ion  battery  accurately. 
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